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Abstract: We present a C++ library for the numerical evaluation of one-loop virtual corrections 
to multi-jet production in massless QCD. The pure gluon primitive amplitudes are evaluated using 
NGLUON [1]. A generalized unitarity reduction algorithm is used to construct arbitrary multiplicity 
fermion-gluon primitive amplitudes. From these basic building blocks the one-loop contribution to 
the squared matrix element, summed over colour and helicities, is calculated. No approximation 
in colour is performed. While the primitive amplitudes are given for arbitrary multiplicities we 
provide the squared matrix elements only for up to 7 external partons allowing the evaluation of the 
five jet cross section at next-to-leading order accuracy. The library has been recently successfully 
applied to four jet production at next-to-leading order in QCD [2]. 
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1 Introduction 

Next-to-leading order (NLO) corrections to multi-particle Standard Model processes at hadron col- 
liders are essential for quantitative predictions of the large QCD background to many new physics 
scenarios. 

Though the computation of the one-loop virtual corrections with four or more final state par- 
ticles has always presented a theoretical challenge, new methods like unitarity [3, 4], generalized 
unitarity [5-8] and integrand reduction [9] have been developed into effective tools suited for a 
wide variety of processes. We also note that the traditional approach — the powerful work horse of 
the last three decades for perturbative calculations in quantum field theory — has been refined and 
improved continuously. For recent developments we refer for example to Ref. [10]. 
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A number of automated approaches have been developed in the last five years [11-16] some 
of which are now publicly available. On-shell techniques have also been used recently to compute 
compact analytic matrix elements within the popular MCFM program [17-22]. The availability of 
one-loop corrections combined with real radiation contributions have made an increasingly large 
number of phenomenological studies possible [23^-6]. More recently improvements beyond fixed- 
order in perturbation theory using parton shower matching techniques [47-49] have been achieved 
for a number of non-trivial cases [50-59]. Despite the tremendous progress there is still room for 
improvement in terms of efficiency and flexibility before one-loop computations reach the same 
level of automation as tree-level predictions and a number of interesting new directions continue to 
appear [10, 60, 61]. 

In this paper we present the C++ library NJet for the numerical evaluation of one-loop am- 
plitudes for multi-jet production in massless QCD. The on-shell generalized unitarity algorithm 
employed in NJet is well suited to deal with the large amount of gluonic radiation which compli- 
cates a traditional Feynman approach to such processes. Our paper is organized as follows: first 
we briefly review the method of generalized unitarity used for the construction of primitive am- 
plitudes with arbitrary numbers of fermion pairs. We then describe the construction of full colour 
and helicity summed squared amplitudes by partial amplitudes in terms of the computed primitive 
amplitudes. A detailed description of the installation procedure and the usage of the NJet package 
is presented in section 3 including a short summary of the interface via the Binoth Les Houches 
Accord (BLHA) [62]. In section 4 we discuss the validation of the NJet library. In particular we 
present cross checks with available analytic results and show a numerical comparison for selected 
benchmark points. In addition we present some performance tests of the program that investigate 
numerical accuracy and speed of both primitive and colour summed amplitudes. 

2 Methods 

The NJet library is built upon the framework of NGLUON [1] which uses generalized unitarity to 
construct primitive gluon amplitudes. This procedure allows one-loop amplitudes to be constructed 
from tree-level amplitudes via the solution to a system of cut equations. One important aspect of 
this method is that it maps the problem of computing any one-loop amplitude to the solution of a 
linear system of equations which is well suited for numerical evaluation. Furthermore the method 
uses on-shell amplitudes instead of individual Feynman diagrams avoiding the significant growth 
in complexity when one-loop amplitudes for large multiplicities are considered. 

The implementation in NGLUON uses the algorithm of Ellis, Giele, Kunszt and Melnikov 
[6, 7] to construct the coefficients of the scalar one-loop integrals. For the rational terms we map 
the problem of cuts in 4 — 2s dimensions to that of four dimensional massive cuts. The scalar one- 
loop integrals are all known analytically and are evaluated using the FF/QCDLoop package [63, 64]. 
An alternative implementation by van Hameren now allows for full complex masses and quadruple 
precision evaluation of the integral functions [65]. 

2.1 Multiple Fermion Primitive Amplitudes 

A new feature added to the NGLUON package to enable the application to multi-jet production 
at NLO is the ability to compute primitive amplitudes with an arbitrary number of external and 
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internal quarks. First results have been presented recently in references [66, 67]. In the following 
we briefly comment on this extension of NGLUON. 

An «-particle, one-loop primitive amplitude is defined by an ordered set of n external particles 
and a set of propagators specifying the loop content. These primitive amplitudes are gauge invariant 
sums of colour ordered Feynman diagrams as introduced in Ref. [68]. The full set of independent 
primitives falls into two subclasses: 

• Amplitudes with a mixture of quark and gluon propagators, A„ , 

• Amplitudes with an internal fermion loop, 

Each primitive amplitude can be associated uniquely with a "parent diagram" [69]. This is a dia- 
gram with n propagators and n colour ordered three-vertices from which all contributing ordered 
Feynman diagrams can be obtained by "pinching" and "pulling" operations. The parent diagrams 
of the two classes of primitives above are determined by the external legs and the parton type of 
the propagator appearing immediately before the first leg. All other internal parton types will then 
follow from the preceding ones, two examples of this are shown in Figure 1. In cases where the 
corresponding vertices do not exist and no internal propagator can be connected we attribute a spe- 
cial propagator as placeholder in order to ensure that the parent diagram always contains a total of 
n propagators. 

Once the parent diagram of the primitive amplitude has been constructed the generalized uni- 
tarity algorithm can be applied to compute the coefficients of the scalar integral basis by factorizing 
each cut into a product of tree-level amplitudes. These amplitudes are then computed with off-shell 
Berends-Giele recursion relations [70]. With a suitable caching system for the off-shell currents, 
this approach scales polynomially in complexity at 0(n 4 ) with n being the number of external legs. 
The power of four arises from the four-gluon vertex, which dominates the recursion for large n, 
independent on the number of quark pairs. The inclusion of fermions is rather a matter of book- 
keeping and no conceptual difficulties arise. 

The rational terms are evaluated by using four-dimensional massive cuts to obtain D-dimensional 
cut information [71, 72]. This corresponds to working in the four-dimensional helicity (FDH) 
scheme [73] which is equivalent to dimensional reduction at one-loop. The coefficients of the D- 
dimensional integral basis [7] are efficiently extracted using a discrete Fourier projection on the 
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Figure 1. Constructing parent diagrams for the primitive amplitudes. In case that no internal propagator 
can be connected since the corresponding vertices do not exist a blank propagator is inserted, represented 
graphically as a shaded blob. 
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internal mass [74, 75]. When only concerned with the rational terms it is possible to replace all 
internal D-dimensional gluons with massive scalars and treat D-dimensional fermions with mas- 
sive fermions. The computation of the rational term is therefore a direct copy of the extraction of 
the integral coefficients but with different tree-level amplitudes which are also straightforward to 
compute with off-shell recursion relations. In addition one can use supersymmetry identities [68] 
to find relations between the rational terms of the two classes of primitives A^ and A^ which 
allows for a more efficient construction of the full colour amplitudes. This has been implemented 
for fermion loops in both the pure gluon and single fermion pair amplitudes: 

£(aL w] (1,...,#0-^(1,...,b)) =0, (2.1) 
^(At ] (l<i,2, I), k q ,{k+l)...,n) 

+ {-\) n A [ :\\ 7l ,n,...,{k+\),k q ,{k-\),...a) 
-AH\\ q ,2,...,(k-\),k q ,(k+\)...,n) 

-(-l) n A\{\l q -,n,...,(k+l),k q ,(k-l),...,2)) =0, (2.2) 
where the operator denotes taking the rational part of the amplitude. 
2.2 Partial Amplitudes and Colour Summation 

To build the full colour amplitudes from the primitive one-loop amplitudes provided by NGLUON 
and its extension to massless quarks we must construct the coefficients of the SU (N c ) colour factors, 
the partial amplitudes. In general the full colour L-loop amplitude can be written as: 

d L \{ Pi }) =£r c (Z ' ) ({a ! })Ag(pi,..., J Pn), (2.3) 

c 

where A„ ;c are the partial amplitudes and the summation is performed over all different colour 
structures T c ({#;}) which depend on the colour of the external partons abbreviated here by a,-. 
For the evaluation of the colour summed squared matrix elements we construct colour matrices C 
by taking the product of colour structures and summing over the colour degrees of freedom of the 
external partons, 

T^a*}), (2.4) 

i 

C^=Z,T c {0 \{a i }).Ty\{a i }). (2.5) 

i 

The colour summed squared Born amplitude will then be constructed from the tree-level partial 
amplitudes A^} according to 

^-Ag^-Ag, (2.6) 
The one-loop interference term is given in complete analogy by 

d^-AW-cg-A®. (2.7) 
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The entries of the colour matrix are given in terms of products of SU (N c ) generators summed over 
colour indices which we evaluate directly into a rational function in N c before hard coding the result 
into C++. 

The algorithm for determining the decomposition of the partial amplitudes in terms of prim- 
itive amplitudes has been constructed along the lines of that suggested by Ellis, Giele, Kunszt, 
Melnikov and Zanderighi [69, 76] and that of Ita and Ozeren [77]. The principle is rather simple: 



1. First we generate all the Feynman diagrams D, for a given «-parton amplitude, separating 
each one into a colour factor denoted by C, and a kinematic factor containing the Lorentz 
structure denoted by Kf. 

■a. = £A=£QZj, (2-8) 

1=1 i=\ 

where A^ia is the total number of diagrams and $dia is the number of diagrams excluding ones 
with four point interactions 1 . Each diagram's colour factor C, is a known linear combination 
of trace factors T c ({a}) with coefficients F c f. 

Ci = Y,T c {{a})F ci . (2.9) 

c 

For clarity we omit the loop superscript on the colour dressed amplitude J^, for the remainder 
of this section. 



2. Secondly we write all possible primitive amplitudes Pi as combinations of colour-stripped 
diagrams: 

/>• Y,M u K h i€{l,2,...,JVp ri } ) N^Njg+NW (2.10) 

7=1 

Y\n q (n-\)\, n q = 2,4,-.. 

where n q is the number of external quarks. Obviously the set {P,} is overcomplete and /V pr j 
is only an upper bound for the number of primitives. The coefficients Mjj have entries of 
±1 and depending on how the parent diagram representing the ordered primitive ampli- 
tude Pj matches onto the kinematic topology of Kj (1 - even-ordered, —1 - odd-ordered, 
- no match). Figure 2 depicts an example of the matching for the primitive amplitude 
^5^(l?'2 u ,3u,4g,5d), the notion of odd and even is taken with respect to the colour ordered 
Feynman rules (see for instance Ref. [78]). 

If some diagrams are known to have additional non-trivial relations, this information can be 
included at this step to reduce the number of columns in My. For instance Furry 's theorem 
can be used on fermionic loops to obtain additional simplifications. 



The kinematic parts K, are understood to contain contributions from both 3-gluon and 4-gluon vertices. 
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Figure 2. Matching Feynman diagrams to the primitive amplitudes Al ' (1j,2 m , 3^,4^,5,/) and 
(l-j,4g,3u,2 u ,5,j). Each vertex is either ordered (red) or unordered (blue) with respect to the primi- 
tive in question. A diagram contributes with a coefficient — 1 when it contains an odd number of unordered 
vertices and + 1 otherwise. Any diagram which cannot be obtained from the original primitive by "pinching" 
and "pulling" operations will have a coefficient. 

3. The rank of the matrix obtained at the previous step defines the number of independent 
primitive amplitudes, which is usually much smaller than the initial upper bound. 

A> pri = rankM, N pii <N pri , A> pri < A> dia (2.12) 

To solve the equation (2.10) it is convenient to consider the augmented matrix M = [M\— 1], 
where 1 is the Apn -dimensional identity matrix. By putting it into reduced row echelon 
form (RREF), we get a solution of (2.10) as the upper Apn rows of RREF(M). The lower 
Apri — Apn rows will contain the left null space of the matrix M, which encodes relations 
between primitive amplitudes. A solution with a minimal basis of primitives is guaranteed 
by the properties of the RREF operation. 

Npn 

Ki=t B U p h {Pjh pA ^{Pj} Npii - (2-13) 

7=1 

4. After substituting solution (2.13) into (2.8) and expressing the colour factors Q using the 
trace basis T c we can form the partial colour amplitudes in terms of independent primitive 
amplitudes Pj as in (2.3), 

#pri /Vdia 

An = £ T c ({a})A n , c = £ T c ({a}) £ Q cj P h Q cj = £ F ci B u . (2.14) 

c c j=\ i=\ 

We are now free to use the result of the primitive decomposition, eq. (2. 14), to evaluate the colour 
sums in eqs. (2.6) and (2.7) in terms of primitive amplitudes. 
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Though this procedure involves generation of Feynman diagrams, only topologies of graphs 
are required making the application even to a relatively high multiplicity straightforward. Never- 
theless the generic algorithm outlined above can be further improved. The four gluon vertex carries 
no additional information as far as colour is concerned and can be ignored leading to a substantial 
speed up. In general, since we are only looking for one specific solution of (2.10), any subset of di- 
agrams which results in the matrix with the same rank as My will be sufficient to obtain the correct 
form of (2.14). 

The number of independent primitive amplitudes for the pure gluon and single fermion pair 
channels follows the well known structure outlined in references [3, 68]. Our approach to the 
fermion loops exploiting Furry type symmetries differs from that of Ita and Ozeren [77] explaining 
the smaller number of primitive amplitudes we obtain in the multi-fermion channels. Table 1 lists 
the number of primitives for all channels required for < 5 jet production. As an example we show 
the fermion loop contribution to the six-quark amplitude 

A 6g (h,2 u ,h-Ad, 5 f , 6,) = A™ xed + B/ Ag oop 
expressed in terms of 4 primitive amplitudes: 

KT P = ^{ T3 + Ta ~ TlNc ~ ~k T6 ) x aW ( 123456 ) + jj- ( r 5 - 23) x AW(124365) 

+ A- (2T 5 - T 2 N C - ^r 6 ) x ^^(125634) + ^-(T 5 - T a ) x A^(126534) 

where 7} are the basis of colour factors 

T[ = 616832854, T 2 = 814836852, T 3 = 812836854, 

t 4 = 814832856, t 5 = 816834852, r 6 = 812834856. 
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24 


60 


60 


&i(g,g,g,g,g,g,g) 


120 


360 


360 


3k(d,d,g,g,g,g) 


24 


120 


33 


&i(d,d,g,g,g,g,g) 


120 


720 


230 


^k(d,d,u,u,g,g) 


12 


80 


13 


Jl 7 (d,d,u,u,g,g,g) 


60 


480 


75 


J?k(d,d,u,u,s,s) 


4 


32 


4 


J%j(d,d,u,u, s,s,g) 


20 


192 


20 



Table 1 . The number of contributing independent primitive amplitudes at tree-level (A [°1 ) and at one-loop for 
the mixed (AM) and fermion loop (A^) cases. Like flavour amplitudes for multiple fermions are obtained 
by symmetrization and therefore contain larger bases of primitives. 

In cases where the final state has a Bose symmetry (i.e. for gluons in the final state) the 
colour summation can be further simplified. Since the phase space integration is symmetric under 
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permutations of the outgoing gluons, it is possible to work with a 'de-symmetrized' version of the 
squared matrix elements and achieve the symmetrization through the phase space integration [79], 

(n-2)!(/j-l)!/2 
J c =l c'=l 

= ^1 /jLIPS n ''I' £ A^ f - tf )4sym -A^ m (2.15) 

Z ^ c=l c'=l 

where <iLIPS n denotes the Lorentz invariant phase space measure of the n parton final state. In 
the case of the gluonic channels the saving is significant and brings the number of independent 
primitive amplitudes from \ {n — 1)! down to (n — 1). The de-symmetrized colour matrices C^ x ' ,As,ym 
have been computed and implemented in NJet for all gluonic channels, the fermionic channels 
however do not improve to the same extent and are implemented in the standard fashion. 

We note that the number of independent primitives is significantly reduced in the case of the 
leading colour approximation where additional symmetries can be exploited. A complete investi- 
gation of the validity of such an approximation remains for a future study. 



2.3 Cache system for tree amplitudes 

The various primitives contributing to a full-colour helicity summed squared matrix element differ 
with respect to the external helicities, the permutations of the external legs and the loop content. 
On the other hand, every primitive is constructed by sewing together partial tree-level amplitudes. 
Therefore partial trees that agree in momenta, flavour and helicities of both external legs and inter- 
nal legs can be reused for different primitive amplitudes. 

An example is illustrated in Figure 3 for a triangle cut of a six-point amplitude with one exter- 
nal quark line and mixed loop content. Flipping the helicity of the fourth external gluon and permut- 
ing the fifth and sixth external leg leaves the five point partial tree amplitudes A (t^ , 1 ~ , 2+ , 3 g , 1%^) 
unchanged. 




1~ 



Figure 3. Triple cut for two primitives with different helicity and ordering of external legs. The partial trees 
A(£* g , l~,2+,3~,.£^) are the same for both primitives. 
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The relation between external permutation and loop momentum requires some comments. 
Let X € {2,3,4,5} denote the topology of the cut, say bubbles, triangles, boxes and pentagons, and 
9Q = {k\ , . . . , &x_ i } a set of x — 1 external momentum sums which determine the kinematic invariants 
of the cut. Since the on-shell loop momentum for a distinct topology depends on x and on the full 
set A necessary condition to reuse loop momenta, is to find equal sets of external momentum 
sums In the above example, % is the same for both permutations {{1,2,3}, {4}, {5, 6}} = 
{{1,2,3}, {4}, {6,5}}. Analogously all six permutations of the first three legs would again lead to 
the same loop momenta. While the agreement of the external momenta and the flavour and helicity 
structure of all legs is required directly by the partial tree amplitude itself, the loop momentum 
matching depends on the topology and the permutation of all external legs. Hence, the knowledge 
of the cuts, permutations, flavours and helicities allows us to assign every partial tree amplitude a 
unique number which we use in a binary tree to cache those trees efficiently. 

Since one has to sum over the helicities of the partons circulating in the loop to determine a 
specific cut, we store always a block of four amplitudes belonging to the four helicity combinations 
(+, +), (+, — ), (— , +) and (— , — ) of the two partons connecting the amplitude to the loop. 



2.4 Universal Pole Structure and Scheme Dependence 

All basic quantities in NGLUON are computed in the four-dimensional helicity (FDH) scheme. 
For consistency checks and the conversion to the conventional dimensional regularisation (CDR) 
scheme the universal structure of infrared (IR) and ultraviolet (UV) poles is required. The universal 
pole structure is known and can be derived from Refs. [80, 81]. For the colour summed squared 
amplitudes in massless QCD the result in the FDH scheme reads: 

+ finite terms, (2.16) 

where denotes the unrenormalized one-loop amplitude in the FDH scheme, n q is the 

number of external quarks and n g is the number of external gluons and g s is the QCD coupling 
constant. The first coefficient of the QCD beta function is given by Po = 11Afc Nf . The quantity 
| 2 is the colour correlated Born amplitude squared defined by 

\^f\ 2 = (A^\T i T j \A { ®) (2.17) 

where T g = f abc and T q = Tfi and we followed the notation used in Ref. [80] to denote by | J?(°)) the 
leading order amplitude which can be understood as a vector in colour space. In order to provide 
results in the conventional dimensional regularisation (CDR) scheme we apply the conversion as 
determined by the universal IR and UV pole structure [80, 81]. The scheme conversion to CDR 
and the renormalization are both proportional to the Born amplitude. For convenience we list the 
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explicit relations between different prescriptions: 



Re[^.(4^ FDH -4^ FDH )) 



(n g + n q -2) 
2 



Po N c 



(2.18) 



gimp 



e 3 



{n g +n q -2)N c 



n, 



(2.19) 



g 2 1^(0) 1 2 



2 3 



6 



4iV c 




g2|^(0)|2 



(2.20) 



where the subscript 7? denotes renormalized amplitudes. 
3 Installation and Usage 

The code can be obtained from https : / /bitbucket . org/n jet/n jet/ as a tar archive. Alterna- 
tively the source can be obtained from the public git repository using 

git clone https://bitbucket.org/njet/njet.git 

Compilation for a wide variety of architectures and compilers is supported via a configure script 
generated with autoconf. The default set-up will be configured by simply using 

. / configure 

which will scan your system for the required pre-requisites. Installation into a local directory 
structure is supported through the option — pref ix=<dir>. Further useful options supported by 
the configure script are: 

• — enable-demo Will compile some additional example programs which can be run with 
make check after the library has been successfully compiled. These demo programs illus- 
trate the usage of the program and the correctness of the implementation. A sample output 
is given in appendix A. 

• — enable-5 jet Includes in the library the one-loop corrections for five jet production. Since 
the compilation requires a significant amount of memory the one-loop corrections for five jet 
production are disabled by default. It is not recommended to use this option with less then 
4 GB of RAM. 

• — disable-cache Per default the cache system is enabled. It allows caching of partial tree 
amplitudes to optimise the colour and helicity sums as described in section 2.3. Due to the 
way the cache is implemented an architecture dependent restriction is imposed on the total 
number of external particles. This is 6 on a 32-bit system and 13 on a 64-bit system. By 
using — disable-cache the cache can be disabled and the restrictions can be lifted however 
at the cost of losing some speed. 
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A complete list of possible options can be obtained using — help. Further details can be found in 
the INSTALL file in the root directory of the distribution. Once the configure script has successfully 
generated the makefiles compilation and installation can be performed with: 

make 

make install 
Additionally, 

make check 

will run some basic tests to verify a correct installation. NJet is distributed with a version of 
QCDLoop for the evaluation of the scalar one-loop integrals and qd [82] to extend the precision of 
the floating-point arithmetic to quadruple precision. 

3.1 Interface via the Binoth Les Houches Accord 

Though the library can be interfaced directly, the simplest way to use it is to link via the Binoth Les 
Houches Accord (BLHA) interface [62]. Once correctly installed, a python script njet.py will 
interpret a user supplied "order file" to produce a "contract file" allowing the various processes and 
options to be passed to NJet. Several sample programs demonstrating the usage of this interface 
are included in the distribution with corresponding order files in the sub-directory examples. These 
examples can be compiled independently using the provided sample Makefile. In the following we 
briefly illustrate an example. The order file used in connection with the BLHA interface consists 
of two parts. The first part contains general settings such as matrix element type and regularisa- 
tion scheme. The second part is a list of elementary sub-processes specified using PDG particle 
codes [83]. An example file for uu — > 3-jets is shown below: 



# order file for u ubar 


-> jets 


MatrixElementSquareType 


CHsummed 


CorrectionType 


QCD 


IRregularisation 


CDR 


AlphasPower 


3 


# process list 




1 -1 -> 21 1 -1 




1 -1 -> 21 2 -2 




1 -1 -> 21 21 21 





Assuming the order file is saved with the file name OLE_order . lh the command 

njet.py -o OLE_contract . lh OLE_order.lh 

will process it and creates a contract file OLE_contract . lh. The user should always check the 
generated contract file to ensure that all requested parameters and sub-processes have been recog- 
nized by NJet which is indicated in the generated contract file. Having a valid contract file one 
can use it with the standard BLHA interface functions which are declared in the n jet . h header: 
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void OLP_Start (const char* filename, int* status); 
void OLP_EvalSubProcess (int men, double* pp, double mur, 

double* alphas, double* rval) ; 

The conventions of their usage are described in [62], we outline here only the main features. Firstly 
the interface has to be initialized by calling the function OLP_Start, passing the file name of a 
contract file as a parameter. After the call, the second parameter status will contain 1 indicating 
the success or indicating the failure of the initialization. In case of a successful initialization 
one can begin using OLP_EvalSubProcess to evaluate the subprocess men for a given phase space 
point with momenta pp, renormalization scale mur and strong coupling alphas. The integer men 
is the number of the process in the contract file (stalling from 1). The momenta are passed through 
the one dimensional array pp. Note that in addition to the momentum for each partem also its mass 
is transferred, e.g. 

pp[5n] = {po,po,pl,pl,m ,p l ,p\,p},p],m l ,...,p n _ l ,pl_ l ,p 2 n _ l ,pl_ l ,m n - l }, (3.1) 

where p' n denotes the i-th momentum component of the n-th parton. The result is returned through 
the array rval with the following order of elements: l/s 2 -pole, 1 /s-pole, finite part, tree amplitude. 
If requested in the order file, the rval array will contain three additional elements: relative accuracy 
of the 1 /s 2 -pole, relative accuracy of the 1/s-pole, and relative accuracy of the finite part. 

3.1.1 Additional BLHA options 

In addition to the standard options (see Ref. [62]) njet .py allows the following additional options 
to control the interface. 

• NJetType [default = loop] 

loop, loopds or tree may be specified to determine the type of the squared matrix element 
returned, loop and tree are conventional colour summed squared amplitudes for one-loop 
and tree contributions respectively, loopds returns the de-symmetrized sums if available 
otherwise the standard one-loop sum is returned. A more advanced usage is to specify a list 
separated by a '-' of a process number followed by a '.' and the type required for that process 
(see example in supplied demo program examples/DSYMtest). 

• NJetReturnAccuracy [default = no] 

If set to yes, the scaling test is used to estimate the accuracy and the result is returned in the 
output array making it three entries longer. 

• NJetSwitchAcc [default = le-5] 

Used to set the relative accuracy at which higher precision floating point arithmetics will 
be used to evaluate the phase space point if double precision has not produced the desired 
accuracy. To account for the statistical nature of the scaling test it is recommended to specify 
1 or 2 digits more than the required accuracy. 

• NJetNf [default = 5] 

Set the number of light quark flavours circulating in the fermion loops. 

A full documentation of the usage of njet .py can be found in lh/LH.txt. 
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3.2 The NParton library for multi-parton primitives 

Also provided with the package is the NParton class for the evaluation of arbitrary multi-parton 
primitive amplitudes both at tree-level and at one loop. The usage is demonstrated in an example 
code compiled into the src-nparton directory: NParton-test . cpp. 

3.2.1 Initialization and Namespaces 

NParton is designed such that both double and quadruple precision can be used in parallel in 
the same program. In order to achieve this, different namespaces NJet_sd and NJet_dd are 
used for single double and double double (^quadruple) respectively. Inclusion of the header file 
src-nparton/NParton-multiprec . h will take care of the necessary declarations. 

Both classes require loop integrals and other parameters to be initialized before they can be 
used which is done by creating the static objects 

static NJet_sd: : Initialize global_init_sd; 
static NJet_dd: : Initialize global_init_dd; 

before the start of the main program. 

3.2.2 Constructor 

The calculation of the primitive amplitudes is hidden in the class NParton. The class constructor 
only requires one integer argument specifying the number of external legs. The namespace for the 
required accuracy must also be specified, e.g. 

NJet_sd: : NParton primitive (int npartons) ; 

for double precision and, 

NJet_dd: : NParton primitive (int npartons); 
for quadruple precision. 

3.2.3 Member functions 

The actual calculation is controlled through member functions of the NParton class. The relevant 
functions are: 

• void setAmp (double moms[][4], int[] flavs) 

This function sets the momenta and flavour of the external legs being the same for all he- 
licity configurations of one primitive amplitude. The four-momenta are entered as an ar- 
ray of four- vectors with the energy in the zeroth component according to (E,p x ,p y ,p z ). 
The list of flavours is specified as follows: for a gluon, (— )i for each (anti-)quark with 
i € {1,2, ...,Nqq} and N q q the total number of external quark lines. Due to fermion and 
flavour number conservation the quarks must always appear as quark-anti-quark pairs. Each 
quark line must be labelled with a different integer. The primitive amplitudes calculated thus 
correspond to the case that all quark flavours are different. The equal flavour case can be ob- 
tained from these results through an appropriate anti-symmetrization. Since QCD conserves 
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the flavour number the fermion lines are not allowed to cross in the primitive amplitudes. A 
valid choice would be, for example, {1,0,2,-2,-1} while {1,0,2,-1,-2} would however 
give an error. 

void eval(int[] hels, int ptype) 
void eval(int hels, int ptype) 

After a previous call of setAmp this function evaluates a single primitive amplitude. The first 
argument is a list of helicities specified by ± 1 for each parton. Alternatively the helicity can 
be given as a single integer between and 2" — 1 which is very convenient to loop over. The 
second argument selects either the mixed or the fermion loop primitive amplitude using the 
internally defined constants NParton: : mixed and NParton: :qloop. 

std: : complex<double> evalTree (int [ ] hels) 
std: : complex<double> evalTree(int hels) 

Similar to the eval function except that only the tree amplitude is evaluated. As eval this 
function requires a previous call of setAmp. The value of the tree-level helicity amplitude is 
returned. 

std: : complex<double> getAmp(int p) 
std: : complex<double> getAbsErr (int p) 
std: : complex<double> getRelErr (int p) 

After calling the eval function the various components of the amplitudes and their absolute 
and relative numerical uncertainty can be accessed with these functions. The different con- 
tributions are obtained by setting the variable p to one of the predefined constants: 
NParton : : epsO returns the finite part, 

NParton : : epsl returns the coefficient of the 1/e pole, 

NParton : : eps2 returns the coefficient of the 1/e 2 pole, 

NParton : : tree returns the tree-level amplitude, 

NParton : : rational returns the rational part, 

NParton : : cutconstructible returns the cut-constructible part. 

void IRcheck() 

After calling the eval function this function will check the poles against known universal 
pole structure. The relative error of both 1/e and 1/e poles compared with the analytic 
formulae will be displayed. 

void refine (double moms, dd_real moms_dd) 

This function converts a double precision phase-space point to quadruple precision to enable 
the calculation in quadruple precision using NJet_dd: : NParton. 

void setMuR2 (double MuR2) 

Used to set the renormalization scale squared, n 2 R . 

std::string ProcessStr() 

Returns a string with the current flavour and helicity structure. This gives the state of the last 
evalAmp call. 
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Process Reference 

pp + 2j (all 2—^2 with colour summation) [84] 

5g [85] 

qq + 3g [68] 

n(g) (finite) [86] 

n(g) (MHV) [87] 

qq + n(g) (finite) [88] 



Table 2. Cross checks against analytic expressions known in the literature. 

With the exception of the refine function, when using NJet_dd: :NParton, double should be 
replaced with dd_real. 

4 Results 
4.1 Validation 

To assess the correctness of our implementation we performed a series of highly non-trivial tests. 
The first test is to check the universal structure of the IR and UV poles (see section 2.4). Both at 
the level of primitive amplitudes and colour summed amplitudes we find for all the processes we 
have tested perfect agreement for the 1/e 2 and 1/e poles where D = 4 — 2s with D denoting the 
dimension of space-time. As far as the UV poles are concerned this check tests a linear combina- 
tion of the coefficients of the two-point scalar integrals. The correct IR structure is sensitive to the 
coefficients of all IR divergent scalar three- and four-point integrals. Since the coefficients of the 
scalar integrals are all calculated the same way reproducing the correct IR and UV structure pro- 
vides a strong test on the entire cut-constructible part. 2 However the rational part of the amplitude 
cannot be tested in this way. To verify the non-trivial rational part a large number of cross checks 
against known processes at individual phase-space points has been performed. We have confirmed 
the results for primitive amplitudes for specific helicity configurations with the known analytic re- 
sults listed in table 2. The checks of the primitive amplitudes in pure gauge theory with the all n 
multiplicity formulae in [86, 87] have already been presented in Ref. [1]. Primitive amplitudes with 
one external quark pair have been tested against the explicit formulae for qq + 3g given in Ref. [68]. 
We find complete agreement for all independent primitive amplitudes and all independent helicity 
configurations for these five-point amplitudes. The IR finite amplitudes with a single quark pair, 
An U ^\qi ,g2 , ■■■■,qt , ■■■gi), have been computed analytically in Ref. [88] to all multiplicity. Direct 
comparisons over all possible configurations of the quark pair have been performed and in all cases 
we find good numerical agreement. The program NParton-qqchecks in the directory src-tools 
performs the comparisons with the analytic formulae listed above. In addition — with the exception 
of the six-gluon channel — we have cross checked numerically the colour summed interferences 
for all processes contributing to three and four jet production with results generated with HELAC- 
1LOOP [90]. We were also able to confirm results for the fermion loops with those generated from 
GoSam [16]. These numerical checks against existing codes can be accessed through the sample 

2 We note that a possible small deviation of the numerically computed poles from the predicted pole structure due to 
numerical rounding errors could give some indication of the numerical accuracy of the current evaluation. 



-15- 





NJet 


BlackHat 


gg ->■ gggg 


born 


+4.9769357371794838 x 10 8 


+4.976935736 x 10 8 


l/£ 2 


-1.7999999999999108 x 10 1 


-1.800000000 x 10 1 


1/8 


-6.5144862057710100 x 10 1 


-6.514486205 x 10 1 


finite 


-3.2130366101334992 x 10 1 


-3.213036625 x 10 1 


dd ->■ gggg 


born 


+2.1622011190045194 x 10 5 


+2.162201118 x 10 5 


l/£ 2 


-1.4666666666663076 x 10 1 


-1.466666667 x 10 1 


1/8 


-5.8264471151950865 x 10 1 


-5.826447114 x 10 1 


finite 


-4.3957884552089730 x 10 1 


-4.395788455 x 10 1 


— )• WWgg 


born 


+1.3745823177248822 x 10 4 


+1.374582317 x 10 4 


l/£ 2 


-1.1333333333333261 x 10 1 


-1.133333333 x 10 1 


1/8 


-4.8019618061344993 x 10 1 


-4.801961805 x 10 1 


finite 


-3.7348157184728159 x 10 1 


-3.734815718 x 10 1 


dd —> uuss 


born 


+3.6761219414819656 x 10 1 


+3.676121941 x 10 1 


l/£ 2 


-8.0000000000001386 


-8.000000000 


1/8 


-2.6337813992821804 x 10 1 


-2.633781399 x 10 1 


finite 


-6.7242846898320696 x 10" 1 


-6.724284689 x 10 1 



Table 3. Numerical comparison between NJet and the 2—^4 amplitudes given in Ref. [89]. 

program NJet-demo. Sample output is provided in appendix A. In case of four jet production full 
agreement has also been found with results for a single phase-space point provided by BlackHat 
[89]. The explicit outcome is shown in Table 3. The results in Ref. [89] are given up to ten decimal 
digits. Apart from the gluonic channel we recover all digits using double precision. For the pure 
gluon channel the error estimate of the numerical evaluation obtained via the scaling test gives, 

A (i),finite = _ 3 .2130366101334992 x 10 1 ±3 x 10~ 8 , 

hence all reliable digits are in agreement with the BlackHat result. When evaluated in quadruple 
precision all digits are in agreement: 

A (i),finite = _ 3 .2130366250275191 x 10 1 ±2 x 10~ 13 . 

The full set of quadruple precision checks can be reproduced using NJet-demo-dd — BlackHat4 j. 

It is important to note that the scalar one-loop integrals are evaluated in double precision which 
limits the final accuracy of NJet to 16 decimal digits even when quadruple precision is used. 

Given the large number of checks passed by NJet we are confident the library produces reli- 
able results for one-loop amplitudes in massless QCD. 
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Since the amplitudes for virtual corrections to five-jet production have not been presented in 
the literature before, we give the numerical values for the evaluation at a single phase space point 
in appendix B. 

4.2 Accuracy 

Given the complexity of the calculations performed in the NJet library it is important to have a 
precise understanding of the numerical accuracy. We have performed a number of tests on the out- 
put of both primitive amplitudes and colour summed squared matrix elements. By performing the 
scaling test as described in Ref. [1] we can reliably estimate the number of correct digits returned 
by the program. 

Figure 4 shows for four channels contributing to four jet production the distribution of the 
estimated accuracy. In total 10 4 phase space points were generated. The accuracy is defined by 



where p are the external momenta, (Jr is the renormalization and x is the scaling parameter. This 
definition gives a direct measure of the valid digits. The phase space points were generated with 
a sequential splitting algorithm as described in [91]. For large multiplicities the algorithm prefers 
collinear configuration. Compared to a flat phase space generation we believe that this simulates 
better what will actually happen in a phase space integration. A weak phase space cut of \pt ■ pj\ > 
s x 10~ 2 was made for a centre-of-mass energy of s = 7 TeV. For the most complex channels around 
1% of points must be re-evaluated in quadruple precision. Some of them are related to small Gram 
determinants but these are not the only sources of numerical instability within the algorithm. As 
we consider channels with more quarks, and therefore lower maximum tensor rank, the accuracy 
improves considerably. For these processes all evaluations pass the threshold accuracy of 10~ 4 after 
quadruple precision has been used. Similar, though significantly improved behaviour, is observed 
for the lower multiplicity channels, plots of the distributions can be found in Ref. [67]. 

Figure 5 shows an identical analysis for two channels contributing to five jet production. For 
the seven gluon amplitude we observe an increasing number of instabilities such that the time to 
re-evaluate with quadruple precision becomes noticeable. Again increasing the number of fermion 
pairs improves the situation a lot. Performing the floating point arithmetics in quadruple precision 
99.99% of points achieve the required 4 digits of accuracy. It might be required switching to octuple 
precision when performing phase-space integration of virtual amplitudes of this complexity. 

4.3 Speed Tests 

An advantage of the recursive construction of ordered tree-level amplitudes within NGLUON is that 
the computation of multi-fermion primitives scales polynomially with time. The scaling properties 
of the algorithm for up to three quark pairs with fixed partem ordering and up to 13 external legs, the 
maximally allowed number with the cache system turned on, are shown in Figure 6. For simplicity 
we restrict the analysis to amplitudes with neighbouring quark pairs, A n {q\,q u ...,qk,qj l ,g2k+\,--;g n ) 
for k quark lines. The times are averaged over all 2" helicity configurations including amplitudes 
that vanish trivially by helicity conservation along fermion lines. The scaling behaviour reflects that 
the cache system has been used to improve the efficiency of the helicity sum and varies between 




(4.1) 
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Figure 4. Accuracy for 4-jet amplitudes: (a) shows the six gluon process, (b) the dd — > 4g process, (c) 
dd — > uu + 2g and (d) the six unlike quark pair process. The thicker histograms show computations in 
double precision whereas the thinner curves show the distribution in quadruple precision for points which 
did not pass the relative accuracy of 10~ 4 when calculated in double precision (region marked by the shaded 
area). Red histograms show the \ poles, green histograms the i and blue histograms the finite part of the 
amplitudes. 



n 4 5 and n 6 . In the asymptotic region the algorithm employed in NGLUON is known to scale as n 8 
for the mixed loop content [1]. For the closed quark loop this is reduced to n 1 since the coupling 
of external off-shell currents to the fermion loop goes exclusively via three-point vertices. We find 
that the asymptotic scaling behaviour is reached at extremely high multiplicities (> 20 legs) and 
has little practical relevance. The asymptotic scaling holds also for the average of configurations 
with non-neighbouring quark pairs. Of course, the timings of individual amplitudes are sensitive to 
a non- vanishing quark antiquark separation. One observes the expected pattern of a speed up in the 
closed quark loop case and a slow down in the mixed quark gluon loop. A more detailed analysis 
for non-neighbouring quark configurations is given in [66]. 

Table 4 shows the timing estimates for colour and helicity summed virtual corrections for all 
sub-processes of multi-jet production with up to 5 jets. The times are quoted for single double 
precision, T se i, and including corrections using quadruple precision in order to achieve an overall 
accuracy of 10~ 4 , 74 digits - We also show the fraction of points which did not pass the relative 
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Figure 5. Accuracy for 5-jet amplitudes: (a) shows the seven gluon process and (b) the dd — >■ dd + 3g pro- 
cess. The thicker histograms show computations in double precision whereas the thinner curves show the 
distribution in quadruple precision for points which did not pass the relative accuracy of 10~ 4 when calcu- 
lated in double precision. Red histograms show the \ poles, green histograms the i and blue histograms 
the finite part of the amplitudes. 




4 5 6 7 8 9 10 11 12 13 

Number of partons 



Figure 6. Estimated evaluation time for primitive amplitude with 0, 1, 2 and 4 fermion lines averaged over 
the total number of helicities as a function of the number of external legs. Both mixed quark gluon loop and 
closed quark loop are shown. The curves are a fit to the polynomial /(«) = an b . The exponent b can be read 
off from the legend. The times were measured on an Intel(R) Core(TM)2 Duo CPU E8400 @ 3.00GHz. 

accuracy when evaluated in double precision. All times include the two evaluations necessary 
to obtain the accuracy estimate via the scaling test. Times have been obtained using the clang 
v2.8 compiler. The time estimated for the reevaluation using quadruple precision is based on the 
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process 


T ,M 
1 sdL i >l 


■* 4 digits l>J 


/V 11ACU J 


UlULCaa 


1 sd l>J 


^4 digits L s J 


(% fixed) 






4g 


u.uju 


U.UJ>U 


\\J.\JVJ ) 




n 99 

u.zz 


99 
u.zz 


(0.22) 






2u2g 


ft frt9 
U.UjZ 


n n^9 

U.UjZ 


^U.UU ) 


ZUJg 


U.J^ 


U.J J 


(0.06) 






2u2d 


n ni 1 

U.Ul 1 


n ni 1 

U.Ul 1 


^U.UU ) 


9n9H 1 rr 
ZUZQlg 


n 1 1 

U. 1 1 


n 1 1 

U. 1 1 


(0.00) 






4u 


O 099 


O 099 


^U.UU ) 


4Ulg 


n 99 
u.zz 


n 99 
u.zz 


(0.03) 




process 




?4 digits [ s ] 


(% fixed) 


process 


7^[s] 


Ta digits [s] 


(% fixed) 


6g 


6.19 


6.81 


(1.37) 


7g 


171.3 


276.7 


(8.63) 


2u4g 


Tin 

/.iy 


H A A 
/.4U 


(U.Jo) 


2u5g 


1 nc 1 


1/11 o 


(3.25) 


2u2d2g 


2.05 


2.06 


(0.08) 


2u2d3g 


45.7 


48.8 


(0.88) 


4u2g 


4.08 


4.15 


(0.21) 


4u3g 


92.5 


101.5 


(1.29) 


2u2d2s 


0.38 


0.38 


(0.00) 


2u2d2slg 


7.9 


8.1 


(0.23) 


2u4d 


0.74 


0.74 


(0.00) 


2u4dlg 


15.8 


16.2 


(0.29) 


6u 


2.16 


2.17 


(0.02) 


6ulg 


47.1 


48.6 


(0.41) 



Table 4. Timing estimates in seconds for full colour and helicity summed virtual corrections. T sl i is the 
time for evaluation in double precision, Z4 digits is me average time estimated to obtain a result correct to 4 
digits using the phase space cuts of section 4.2. All times include the two evaluations necessary to obtain the 
accuracy estimate via the scaling test and were obtained on an Intel(R) Xeon(R) CPU E3-1240 @ 3.30GHz. 







gg^lg 


gg^4g 


gg^5g 


standard sum 


0.03 


0.22 


6.19 


171.31 


de-symmetrized 


0.03 


0.07 


0.57 


3.07 



Table 5. Timing estimates in seconds for the de-symmetrized colour and helicity summed gluonic channels. 
All times include the two evaluations necessary to obtain the accuracy estimate via the scaling test and were 
obtained on an Intel(R) Xeon(R) CPU E3-1240 @ 3.30GHz. 

percentage of rescued points from the accuracy tests in the previous section. On average quadruple 
precision evaluation was found to be between 7 and 8 times longer than double precision. We stress 
that even in the worst case of the seven gluon amplitude the required reevaluation only doubles the 
total run time. 

Table 5 shows the evaluation times using the de-symmetrized colour sums for the pure gluonic 
channels which exploit the Bose symmetry of the final state. Though these channels benefit the 
most from this treatment the complex channels with a single fermion pair would also see a con- 
siderable speed up. Again we stress that these de-symmetrized sums contain the same full colour 
information as the standard ones after the integration over the phase-space or equivalently after 
explicitly summing the ^(n — 2)! permutations of the final state gluons. 

5 Conclusions 

In this paper we have presented the C++ library NJet for the numerical evaluation of one-loop 
amplitudes for multi-jet production at hadron colliders. Using generalized unitarity together with 
off-shell recursion relations we were able to construct multi-fermion primitive amplitudes at a 
computational cost growing polynomially with time. Accuracy estimates obtained by applying a 



-20- 



momentum re-scaling test are used to ascertain when higher precision numerics are required to 
achieve a preset number of correct digits. A simple cache system allows the efficient computation 
of helicity and permutation sums. 

In order to use the multi-fermion primitive amplitudes for NLO QCD computations the colour 
space information must be included. By matching primitive and partial amplitudes via Feynman 
diagrams we obtained full colour formulae for all channels contributing to < 5 jet production. We 
employed Furry type symmetries and linear algebra techniques to minimize the number of primitive 
amplitudes appearing in the colour sums. Furthermore we made use of Bose symmetries in the final 
state to simplify the most complex all-gluon processes. Since the on-shell algorithm is particularly 
well suited to the computation of multi-gluon states which contain degenerate degrees of freedom 
we find a favourable performance in comparison to existing public one-loop amplitude generators. 

The interface to the library is either directly to the primitive amplitudes using the NParton 
class or via a Binoth Les Houches Accord [62] to the colour and helicity summed amplitudes. The 
latter allows the library to be linked with existing parton level Monte Carlo programs which can 
combine the virtual and real radiation contributions to obtain full NLO predictions. The package 
presented here has been recently used to calculate three- and four jet cross sections at next-to- 
leading order in QCD [2]. This application nicely illustrates that the tool presented can be used for 
state of the art calculations. We also provide in Table 6 in the appendix benchmark points for all 
one-loop corrections required for the evaluation of the five jet cross section at next-to-leading order 
accuracy. 

Future developments to include heavy quark corrections would be interesting. 

The library NJet is publicly available from the project home page at https : / /bitbucket . org/n jet/n jet . 
Acknowledgments 
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A Numerical cross checks 

With — enable-demo used to configure the package two example programs are compiled into their 
respective source directories: 

src-tools/N Jet-demo 
src-nparton/NParton-demo 

NParton-demo is a short program showing how to use the NParton library by evaluating some one- 
loop multi-parton primitive amplitudes. The code found in src-nparton/NParton-demo . cpp has 
some comments to explain the functions calls. 

A number of hard coded cross checks for the colour summed amplitudes can be found in the 
program NJet-demo. The individual checks can be accessed with command line options, simply 
type NJet-demo — help to see the available options. 
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For convenience we list the program output for the comparison with HELAC-1LOOP. 

NJet-demo — Helacll gives a numerical cross check with HELAC-1LOOP [90] for the 
available channels of dijet, three-jet and four-jet production. We find again good numerical agree- 
ment. The exact phase-space points for these checks are displayed to double precision when run- 
ning the program. The results of these checks are as follows: 



==== — > g g g g ==== 
NJet: 

born: +1 . 4064482826874553e+02, 

l/eps"2: -2. 7 96310037 6599282e+01, 
1/eps: -2.4893580169045921e+00, 
finite: +2 . 2538771232905535e+01, 

==== — > u~ u g g ==== 
NJet: 

born: +3 . 8973755683505198e+00, 

1/eps "2: -5.5963473237272909e-01, 
1/eps: -2.2701782897620162e-01, 
finite: +1 . 58 12 14 9573 18 63 10e-0 1 , 

==== — > u~ u d~ d ==== 
NJet: 

born: +5 . 5757038821628091e-01, 

1/eps ~2: -4. 926956607 1365262e-02, 
1/eps: -4.9621618996857207e-02, 
finite: -2 . 8 7 324 24 65 37 63 1 96e- 02 , 

==== — > u~ u u~ u ==== 
NJet: 

born: +5 . 4442906439318520e+01, 

1/eps "2: -4.8108336321631899e+00, 
1/eps: +3.8870984690869954e+00, 
finite: +2 . 3627403557985861e+01, 

2 — > 3 

==== — > g g g g g ==== 
NJet: 

born: +5 . 0536873228500680e-02, 

1/eps "2: -1. 255971936724 6716e-02, 
1/eps: -6.5137043664892195e-03, 
finite: +1 . 9628928574388416e-02, 



Helaclloop : 


rel 


.cliff 


+1.4064482826874587e+02, 


+2. 


4e-15 


-2.7963100376599272e+01, 


+ 3. 


8e-16 


-2.4893580169051024e+00, 


+2. 


Oe-13 


+2.2538771232905539e+01, 


+ 1. 


6e-16 


Helaclloop : 


rel 


.cliff 


+3.8973755683505242e+00, 


+ 1. 


le-15 


-5.5963473237273020e-01, 


+2. 


Oe-15 


-2.2701782897620176e-01, 


+ 6. 


le-16 


+1.5812149573186496e-01, 


+ 1. 


2e-14 


Helaclloop : 


rel 


.cliff 


+5.5757038821628113e-01, 


+ 4. 


Oe-16 


-4.9269566071365241e-02, 


+ 4. 


2e-16 


-4.9621618996857221e-02, 


+2. 


8e-16 


-2.8732424653763248e-02, 


+ 1. 


8e-15 


Helaclloop : 


rel 


.cliff 


+5.4442906439318527e+01, 


+ 1. 


3e-16 


-4.8108336321631695e+00, 


+ 4. 


2e-15 


+3.8870984690871220e+00, 


+ 3. 


3e-14 


+2.3627403557985797e+01, 


+2. 


7e-15 


Helaclloop : 


rel 


.cliff 


+5.0536873228500589e-02, 


+ 1. 


8e-15 


-1.2559719367246622e-02, 


+ 7. 


5e-15 


-6.5137043664457188e-03, 


+ 6. 


le-12 


+1.9628928574411668e-02, 


+ 1. 


2e-12 
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-> u' 


u g g g ==== 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 5 


,0328282104265365e- 


04, 


+5.0328282104265484e- 


04, 


+2. 


4e- 


15 


l/eps~2 : 


-9 


,7283503961992726e- 


05, 


-9.7283503961991818e- 


■05, 


+ 9. 


3e- 


15 


1/eps : 


-6 


,1775045952609201e- 


05, 


-6.1775045952603048e- 


■05, 


+ 1. 


Oe- 


13 


finite : 


+ 1 


,1409871173998674e- 


04, 


+1.1409871173998645e- 


04, 


+2. 


5e- 


15 


— - 


-> u' 


u d" d g ==== 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 3 


,6329475567709317e- 


04, 


+3.6329475567709729e- 


■04, 


+ 1. 


le- 


14 


l/eps~2 : 


-5 


,0160076608990277e- 


05, 


-5.0160076608990379e- 


■05, 


+2. 


Oe- 


15 


1/eps : 


-4 


,1238301409751345e- 


05, 


-4.1238301409752870e- 


■05, 


+ 3. 


7e- 


14 


finite : 


+ 1 


,6929367019718101e- 


05, 


+1.6929367019719798e- 


05, 


+ 1. 


Oe- 


13 




-> u' 


u u~ u g ==== 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 1 


,6607841328194570e- 


03, 


+1.6607841328194737e- 


03, 


+ 1. 


Oe- 


14 


l/eps~2 : 


-2 


,2930432666983780e- 


04, 


-2.2930432666983907e- 


■04, 


+ 5. 


6e- 


15 


1/eps : 


-1 


,6167755068577280e- 


04, 


-1.6167755068572718e- 


■04, 


+2. 


8e- 


13 


finite : 


+ 7 
- 2 - 


,0356660045273211e- 
__> 4 


04, 


+7.0356660045274610e- 


■04, 


+2. 


Oe- 


14 


==== - 


-> u' 


u g g g g === 


== 














NJet: 




Helaclloop : 




rel 


.diff 


born : 


+2 


,4749786907729620e- 


05, 


+2.4749786907729492e- 


■05, 


+ 5. 


2e- 


15 


l/eps"2 : 


-6 


,0142737839968278e- 


06, 


-6.0142737839968761e- 


06, 


+ 8. 


0e- 


15 


1/eps : 


-1 


,4238548051572628e- 


05, 


-1.4238548051576113e- 


■05, 


+2. 


4e- 


13 


finite : 


-2 


,1818575805881747e- 


06, 


-2.1818575805944554e- 


06, 


+2. 


9e- 


12 


==== - 


-> u' 


u d" d g g === 


== 














NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 3 


,8990501849846650e- 


06, 


+3.8990501849846963e- 


06, 


+ 8. 


Oe- 


15 


l/eps~2 : 


-7 


,3214448592467342e- 


07, 


-7.3214448592467173e- 


07, 


+2. 


3e- 


15 


1/eps : 


-1 


,5000972850719066e- 


06, 


-1.5000972850719574e- 


06, 


+ 3. 


4e- 


14 


finite : 


+ 1 


,1570207960122947e- 


07, 


+1.1570207960076906e- 


07, 


+ 4. 


0e- 


12 




-> u' 


u u~ u g g === 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 1 


,2029972817541046e- 


05, 


+1.2029972817541090e- 


05, 


+ 3. 


7e- 


15 


l/eps~2 : 


-2 


,2589291869350247e- 


06, 


-2.2589291869350408e- 


06, 


+ 7. 


le- 


15 


1/eps : 


-4 


,9249280446345875e- 


06, 


-4.9249280446538397e- 


06, 


+ 3. 


9e- 


12 
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finite : 



+3.6301807851586739e-06, 



+3.6301807851561252e-06, 



+7.0e-13 



==== - 


-> u" 


u d~ d s" s === 


: = 














NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 4, 


6602437519749336e- 


08, 


+4.6602437519749931e- 


08, 


+ 1. 


3e- 


14 


l/eps~2 : 


-6, 


1770188741857187e- 


09, 


-6.1770188741856584e- 


09, 


+ 9. 


8e- 


15 


1/eps : 


-1, 


2079244035518674e- 


08, 


-1.2079244035518879e- 


08, 


+ 1. 


7e- 


14 


finite : 


+ 3, 


0749739659342438e- 


09, 


+3.0749739659321969e- 


09, 


+ 6. 


7e- 


13 


==== - 


-> u' 


u d~ d d~ d === 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 1, 


7752256666856891e- 


08, 


+1.7752256666857116e- 


08, 


+ 1. 


3e- 


14 


l/eps~2 : 


-2, 


3530104931548807e- 


09, 


-2.3530104931548365e- 


09, 


+ 1. 


9e- 


14 


1/eps : 


-4. 


8669803368255748e- 


09, 


-4.8669803368265401e- 


09, 


+2. 


Oe- 


13 


finite : 


+ 1. 


2351513905489450e- 


09, 


+1.2351513905470286e- 


09, 


+ 1. 


6e- 


12 




-> u" 


u u~ u u~ u === 
















NJet: 




Helaclloop : 




rel 


.diff 


born : 


+ 1. 


9800262627817160e- 


06, 


+1.9800262627817295e- 


06, 


+ 6. 


8e- 


15 


l/eps~2 : 


-2. 


6244677848455005e- 


07, 


-2.6244677848455651e- 


07, 


+2. 


5e- 


14 


1/eps : 


-3. 


5478617113255096e- 


07, 


-3.5478617113278930e- 


07, 


+ 6. 


7e- 


13 


finite : 


+ 1. 


2984349868781561e- 


06, 


+1.2984349868773965e- 


06, 


+ 5. 


8e- 


13 
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B Numerical evaluation for 7-parton virtual corrections 



We evaluate all indepenent channels contributing to the virtual corrections to pp - 
following phase space point: 

Pi = {-5.0000000000000000 x 10 2 , 0.0000000000000000 x 10°, 

0.0000000000000000 x 10°, -5.0000000000000000 x 10 2 }, 
p 2 = {-5.0000000000000000 x 10 2 , 0.0000000000000000 x 10°, 

0.0000000000000000 x 10°, 5.0000000000000000 x 10 2 }, 
p 3 = {8.6354068143781365 x 10 1 , -1.5213389320261800 x 10 1 , 

3.7633551294916273 x 10 1 , -7.6218722682185415 x 10 1 }, 
p 4 = {2.8011818180937634 x 10 2 , -8.3126111650582232 x 10 1 , 

-2.6320385675865049 x 10 2 , 4.7749085116026578 x 10 1 }, 
p 5 = {1.2752252956966052 x 10 2 , -9.0449041295993482 x 10 1 , 

-8.3178307703078929 x 10 1 , 3.4093043339258053 x 10 1 }, 
p 6 = {4.1413006837454346 x 10 2 , 2.3214556494593859 x 10 2 , 

3.3275443678081871 x 10 2 , -8.2985751852442578 x 10 1 }, 
p 7 = {9.1875152102638314 x 10 1 , -4.3357022679101075 x 10 1 , 

-2.4005823614005564 x 10 1 , 7.7362346079343363 x 10 1 }. 

The result has been obtained with the number of light quark flavour rif = 5, a renormalization scale 
of pr = 91.188 and a s = 0.118. Table 6 shows the numerical values of the amplitudes. 
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